Rotating electrons in quantum dots: Quantum Hall liquid in the classical limit 
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We solve the problem of a few electrons in a two-dimensional harmonic confinement using quantum 
mechanical exact diagonalization technique, on one hand, and classical mechanics, on the other hand. 
The quantitative agreement between the results of these two calculations suggests that, at low filling 
factors, all the low energy excitations of quantum Hall liquid are classical vibrations of localized 
electrons. The Coriolis force plays a dominant role in determining the classical vibration frequencies. 
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Spectroscopic studies of semiconductor quantum dots 
have shown that in zero magnetic field they behave like 
'artificial atoms' showing electronic shell structure and 
obeying Hund's ruleQ, while in a strong magnetic field 
they have properties of quantum Hall liquidsQ, • The 
integer quantum Hall effect corresponds in a quantum dot 
the formation of so-called maximum density droplet |j] of 
polarized electrons. For higher magnetic fields the re- 
search has mainly been theoretical, showing formation of 
quantum liquids related to the fractional quantum Hall 
effect (QHE)H H Q. Although the quantum dot is fi- 
nite and has a surface, its electronic structure at high 
magnetic fields shows properties of the QHE, which is 
usually studied using geometries with periodic boundary 
conditions, i.e. the surface of a sphere or a torusQ. 

The generic model of a semiconductor quantum dot is 
very simple: Electrons interacting with Coulomb interac- 
tion are confined in a two-dimensional harmonic poten- 
tial. Due to the circular symmetry, the angular momen- 
tum is a good quantum number, and due to the harmonic 
confinement the center-of-mass motion exactly separates 
out EH EH from the internal motion of the particles. The 
wave functions of non-interacting electrons are simple 
(with or without an external magnetic field). Moreover, 
in a strong magnetic field, once the electrons are polar- 
ized, the only effect of the magnetic field is to put the 
electron system in rotation, i.e. increase the total angu- 
lar momentum. For small number of electrons the many- 
particle problem can be solved numerically exactly. It is 
then not surprising that the quantum dot geometry has 
been extensively used for testing different models of the 
fractional quantum Hall liquid ( e.g. the Laughlin wave 
function^] and Jain construction^^)- 

The purpose of this letter is to show that the quan- 
tum mechanical many-particle spectrum of strongly cor- 
related electrons of the QHE system can be quantitatively 
determined by solving the Newton equation of motion. 
To this end we solve, for small number of electrons, the 
excitation spectra using many-particle quantum mechan- 
ics and compare the results to those determined by quan- 
tizing the rotational and vibrational modes calculated in 
a rotational frame for a Wigner molecule using classical 
mechnanics. 



We assume a generic model of a polarized (or spin- 
less) electrons interacting with unscreened Coulomb in- 
teraction in a two-dimensional harmonic potential. The 
Hamiltonian is 
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where N is the number of particles, m the electron mass, 
r = (x, y) a two-dimensional position vector, and ujq the 
oscillation frequency of the confining potential. We solve 
the many-particle problem for a fixed angular momentum 
using the single particle basis of the lowest Landau level 
(LLL): 



Mr, <j>)=A e r e e- r2 / 4 e u 
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where £ is the single-particle angular momentum and An 
a normalization factor. We are interested in solutions for 
large total angular momenta. In this case the restriction 
of the basis in the LLL is a good approximation 13]. Note 
that we do not explicitly include magnetic field in our 
calculations: The only effects the magnetic field would 
have, were to polarize the electron gas due to the Zee- 
man effect and to increase the total angular momentum. 
We present our results in terms of angular momentum. 
The total angular momentum can be related to the filling 
factor of electrons with the relation v — N(N — 1)/2L, 
where N is the number of electrons and L the total an- 
gular momentum. 

We solve the many-particle Schrodinger equation us- 
ing the straightforward CI technique. The many-particle 
states are linear combinations of Slater determinants of 
single particle configurations. The number of configu- 
rations needed increases fast with N and L, but for the 
results shown here a complete basis in the LLL could still 
be used. For determining the Coulomb matrix elements 
we used the method of Ref. Q 

Classical electrons in a 2D harmonic trap form Wigner 
molecules, where the electrons arrange in consecutive 
circles [H The elementary excitations of the clas- 

sical system consist of the center of mass vibrations with 
angular frequency u)q, internal vibrations which can be 
solved from the dynamical matrix, and rigid rotations of 
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the Wigner molecule. Since we are interested in the exci- 
tations of a rotating system (with fixed angular momen- 
tum) we have to solve the internal vibrations in a rotating 
frame^3- Taking into account the Coriolis force turns 
out to be crucial in the case of electrons in a harmonic 
trap (although its effect in real molecules is small). 
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FIG. 1: Classical vibrational frequencies as a function of the 
angular velocity of the rotating Wigner molecule with 2, 3, 
and four electrons. The breathing mode (only mode for two 
electrons) is shown as a dotted line. The solid (dashed) lines 
show other modes for four (three) electrons (in atomic units 
with uuo = 0.5). 

The equilibrium positions of the electrons depend on 
the angular velocity uj r or angular momentum L = Icu r 
of the molecule (I is the moment of inertia I = Y] mrf), 
and they can be solved by minimizing the classical energy 
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The eigenfrequencies of the vibrations can then be solved 
analytically from the equations of motion of the rotating 
frame (by linearizing the equations around the equilib- 
rium positions of electrons). Results for classical vibra- 
tional frequencies for 2, 3 and 4 electrons are shown in 
Fig. 1. The eigenfrequencies are shown as a function of 
the angular velocity of the rotation u> r . We notice that 
all frequencies approach to either or 2a;o. In all cal- 
culations shown in this letter we have used atomic units 
(to = h = e = ciq = 1) and chosen luq = 0.5. We should 
note that in the LLL the results are independent of ujq, 
apart from the energy scale. 

For each particle number the highest mode corresponds 
to the breathing mode where the molecule expands and 
shrinks without changing its shape. In a rotating frame, 
however, each electron moves along an ellipse around its 
equilibrium position. In the case of two electrons this 
mode is naturally the only vibrational mode. For 2, 3 and 
4 particles the breathing mode energy is the same and 



its dependence of the angular velocity is similar Wb m = 
\/3ui + bj r (this is general result for a Wigner molecule 
where the electrons form one single ring). 

In the case of 3 and 4 electrons the non-rotating system 
has degenerate vibrational modes. However, when the 
system is put in rotation the degeneracy will split. In the 
case of three electrons the ground state is an equilateral 
triangle. In the doubly degenerate vibrational mode the 
triangle stays isosceles, the base stretching and shrinking. 
In the non-rotating system we can choose any of the sides 
to be the base, but only two of them linearly indepen- 
dently. However, we can also choose a linear combination 
of the two modes in such a way that the base circulates 
around. In this pseudo-rotation the electrons move in 
circles around their equilibrium positions (as indicated 
in Fig. 4). The degeneracy then follows from the two 
possible direction of the pseudo-rotation. If the molecule 
is put in rotation, the Coriolis force splits this degeneracy. 
Similar analysis can be made for the degenerate mode in 
the case of four electrons. It is interesting to note that 
in the case of sodium trimer, where the Jahn-teller effect 
opens the equilateral triangle, a related pseudo-rotation 
have been observed 251. 
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FIG. 2: Many-particle energy spectrum for three electrons. 
The interaction energy is shown as a function of the angular 
momentum. The black dots are results of the Schrodinger 
equation and the open squares the results from the classical 
model. The center of mass excitations are not shown as points 
but are indicated as dashed lines. The numbers indicate the 
order of the vibrational state (ru). 

Having the classical vibrational frequencies in hand, 
we follow the normal assumptions in quantization of 
phonons0 and get an estimate for the total energy as 



^2fuj k (n k + -) + huj (n + 1), (4) 



where io k are all the vibrational frequencies determined 
in the rotating frame and n k = 0, 1, 2, • • • , and the 
last term corresponds to the center of mass excitations. 
These energies can now be compared with those obtained 
by solving the many-particle Schrodinger equation. How- 
ever, we still have to note that the quantum state for 
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FIG. 3: Many-particle energy spectrum for four electrons. 
The interaction energy is shows as a function of the angular 
momentum. The black dots are results of the Schrodinger 
equation and the open squares the results from the classical 
model. The center of mass excitations are not shown as points 
but are indicated as dashed lines. The numbers indicate the 
vibrational state (ni,Ti2). 



fermions has to be antisymmetric. This means that not 
all combinations of L and nk are allowed (but any center- 
of-mass excitation is always possible). The symmetry 
analysis can be made with help of group theory [l3,|20, 21] 
and it shows, for example, that without any vibrational 
modes only every second angular momentum is allowed 
for two particles, every third for the triangle of three par- 
ticles and every fourth for the square of four particles. 
The lowest energy of any other angular momentum value 
must have either a center of mass excitation or an internal 
vibration in addition of the internal rotation. This fact 
leads to the well known oscillations of the lowest energy 
as a function of the angular momentum^, Ull l22| . 

Figures 2 and 3 show the lowest energies for three and 
four particles obtained using two different ways: Solv- 
ing exactly the many-particle Schrodinger equation and 
by solving the vibrational modes using classical mechan- 
ics and determining the energy from Eq. (2). The 
figures show the 'interaction energy', the energy differ- 
ence between the interacting and non-interacting elec- 
trons, E — (N + L)fajj), where E is the solution of the 
Schr':odinger equation of of Eq. (4). In this presen- 
tation the center-of-mass excitations become horizontal 
(the center-of-mass motion is independent of the interac- 
tions) and they are only indicated as dashed lines. The 
results show that spectra determined by quantizing the 
classical vibrations agree very well with those from the 
full quantum mechanical calculations. The error in the 
total interaction energy is less than 1 %. We note in 
passing that in the case of the frequently studied system 
of two electrons j^, 01 the low energy spectrum only 
consists of rigid rotations and center-of-mass excitations, 
since the only vibrations mode is high in energy. 

In the case of three electrons the purely rotational state 
can occur at angular momenta L = 3n, where n is an in- 
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FIG. 4: The left column shows the classical geometries 
and schematically the low energy vibrational modes (pseudo- 
rotations) for three and four electrons. The two other columns 
show contour plots of the pair correlation functions from 
the quantum mechanical calculation. The reference point is 
shown as a cross. The relation of the pair correlation functions 
to the classical extreme geometries is indicated with solid and 
dashed lines. The numbers in parentheses are (L,ni,ri2) in- 
dicate the angular momentum and the vibrational state. 



teger. For L = 3n + 1 the lowest energy state has a 
center of mass excitation (the point not shown) and for 
L = 3n + 2 the lowest energy state has vibrational mode 
u)\. The excited states for any angular momentum show 
simple systematics consisting of center of mass excita- 
tions and multiplets of oj\, as indicated in Fig. 3. 

In the case of four particles the energy spectrum is 
more complicated due to the two low energy vibrational 
modes. Nevertheless, Fig. 3 shows that also in this 
case the spectrum consists of periodic sequences, but the 
length of the period is now four. 
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In Figures 2 and 3 we have shown the spectra for large 
angular momenta, where the quantization of the classical 
energies give quantitatively accurately the energy differ- 
ences. The angular momentum L — 54 for four particles 
and L = 27 for three particles corresponds to to the fill- 
ing factor v = 1/9 of a quantum Hall liquid, where the 
electrons are expected to be localized to a Wigner crys- 
tal. However, the excitation spectrum, when plotted as 
in Figs. 2 and 3 are qualitatively similar already from an- 
gular momenta corresponding to filling factor v = 1/3. 

Finally, we want to show that the pair correlation func- 
tions determined from the full quantum mechanical cal- 
culation are consistent with the vibrational modes deter- 
mined from classical mechanics. Figure 4 shows schemat- 
ically the classical vibrational modes in the the rotational 
frame compared to the pair correlation functions deter- 
mined for different states shown in the spectra in Figs. 2 
and 3. The pair correlation functions for the purely ro- 
tational states shows clearly the localization of the elec- 
trons, and their shape is in agreement with the semiclas- 
sical model ov Matulis and Anisimovas^J- However, the 
interpretation of the vibrational states is not so simple, 
because the structure of the pair correlation function de- 
pends on the reference point. In Fig. 4 the reference 
electron is chosen to be at the extreme point of the clas- 
sical motion. In this case we can identify the classical 
extreme positions as shown as triangles and quandran- 
gles drawn on top of the pair correlation functions. Note 
that the pair correlation function can capture the classi- 
cal geometry in two different positions. 

In conclusion, we have determined the classical ro- 
tational and vibrational modes of electrons in a two- 
dimensional harmonic trap using a rotational frame. Re- 
sults of the quantized energies are compared with numer- 
ically exact solutions of the many-particle Schrodinger 
equation, and they indicate clearly that when the an- 
gular momentum of the system increases to the region 
which corresponds to the filling factor of 1 /9 of the quan- 
tum Hall liquid, the whole low-energy excitation spec- 
trum can be quantitatively described with the classical 
rotations and vibrations. The pair-correlation functions 
are consistent with the classically determined vibrational 
modes. 
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